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ABSTRACT 

We extend the abundance matching technique (AMT) to infer the satellite-subhalo and central- 
halo mass relations (MRs) of local galaxies, as well as the corresponding satellite conditional mass 
functions. We use the observed galaxy stellar mass function (GSMF) decomposed into centrals and 
satellites and the A-CDM distinct halo and subhalo mass functions as inputs. We explore the effects 
of defining the subhalo mass, m su b, at the time of (sub)halo accretion (m^) versus defining it at the 
time of observation (m°^); and we test the standard assumption that centrals and satellites follow 
the same MRs. We show that this assumption leads to predictions in disagreement with observations, 
specially when mj is used. Instead, we find that when the satellite-subhalo MRs are constrained by 
the satellite GSMF, they are always different from the central-halo MR: the smaller the stellar mass, 
the less massive is the subhalo of satellites as compared to the halo of centrals of the same stellar 
mass. This difference is more dramatic when is used instead of m?^. On average, for stellar 

masses lower than ~ 2 x 10 11 M Q , the dark mass of satellites decreased by 60 — 65% with respect to 
their masses at accretion time. We find that MRs for both definitions of subhalo mass yield satellite 
conditional mass functions (CSMF) in agreement with observations. Also, when these MRs are used 
in a halo occupation model, the predicted two-point correlation functions at different stellar mass bins 
agree with observations. The average stellar-halo MR is close to the MR of central galaxies alone, 
and conceptually this average MR is equivalent to abundance matching the cumulative total GSMF 
to the halo + subhalo mass function (the standard AMT). We show that the use of m°^£ leads to less 
uncertain MRs than mj^£, and discuss some implications of the obtained satellite-subhalo MR. For 
example, we show that the tension between abundance and dynamics of Milky- Way satellites in the 
A-CDM cosmogony gives if the faint-end slope of the GSMF upturns to a value of ~ —1.6. 
Subject headings: galaxies: abundances — galaxies: evolution — galaxies: halos — galaxies: luminos- 
ity function, mass function — galaxies: statistics — cosmology: dark matter. 



1. INTRODUCTION 

In recent years the abundance matching tech- 
nique (AMT) has emerged as a simple yet power- 
ful statistical approach for connecting galaxies to ha- 
los without requiring knowledge of the underlying 
physics (e.g.JVale fc Ostrikerl2004l:|Kravtsov et al.| | 200l 



Conrov. Wcchslcr & Kravtsov 2006; Shankar et al. 2006; 
Baldr vTGlazebrook fe Driver! 120081: iConrov fc Wechsler l 
20091: iDrorv et alJ 120091 iBehroozi Conrov fc Wechslerl 
20101 and references therein). 



Briefly, the AMT assumes a one-to-one monotonic 
relationship between stellar and halo masses which 
can be constrained by matching the cumulative ob- 
served galaxy stellar mass function (GSMF) to the 
theoretical halo plus subhalo cumulative mass func- 
tion. Interestingly enough, this simple approach suc- 
cessfully reproduc es the observed spatial clustering 
of galaxies (e.g., lConrov~W cchsler fe Kravtsov! 120061 : 



iMoster et afll2010( ). The AMT allows to probe the av- 
erage galaxy stellar-halo mass relation, M*(Mh) (here- 
after SHMR), delivering very useful informati on for con- 
strain i ng models of galaxy evolution (e.g.. IGuo et all 



20101: iFhnTani, Avila-Reese fc Rodriguez-Pueblal 120101: 



Avila-Reese et aLll2011li " 

The above has motivated several authors to use the 
AMT extensively. For example, with the advent of 



large galaxy surveys at different redshifts, the AMT 
has been applied for c onstraining the evolution of 
the a v erage SHMR ( e.g.. IConrov. Wechsler fc Kravtsoyl 
2001 IDrorv et all 120091: IConrov fc Wechslerl 120091: 
Moster et alll2010UBehroozi. Conrov fc Wechslerll2010D . 
As a natural extension, these studies have been com- 
bined with predicted average halo mass aggregation 
histories in order to infer average galaxy growth 
histo ri es as a function of mass (|Conrov fc Wechslerl 
12001 IFirmani fc Avila-Reesd l20la~ see for a re- 
view Avila-Reese & Firmani 2011, and references 
therein). By including observational information 
on the gas content of galaxies, the AMT has 
been also used to con strain the baryon mass to 
Mh r e lation of galaxie s (Baldry, Glazc brook fc Driverl 
120081: IRodriguez-Puebla et all 120 111) . Finally, vari- 
ants of the AMT, where instead of mass func- 
tions, circular velocity functions or functions of any 
other galaxy/halo global property are employed, have 
been explored, too (e.g.. IConrov Wechsler fc Kravtsov 



2011 ) 



2006; Bla nton. Geha fc West 2008; Truiillo-Gomez et al 



apucblaOast 
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The AMT has been commonly applied to the to- 
tal (central plus satellite galaxies) GSMF matched 
against the total (distinct plus satellite) halo popula- 
tion. This approach has been criticized, because quite 
different average SHMRs are obtained for different pro- 
posed forms of the satellite stellar-subhalo mass relation 
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(SSMR, m.k.b)) and the central SHMR (M*(M h )Q; 
iNeistein et al]|2011[ ). 

A common (questionable) assumption is that the 
SSMR is identical to the central SHMR. Under this as- 
sumption, it is also common to define subhalo mass at 
the time of accretion rather than at the time of 

observation (mj), when subhalos have lost a signifi- 
cant fraction of mass due to tidal stripping. The use of 
m sub nas been justified because this way is avoided the 
question of subhalo mass loss, and regarding the satel- 
lite to*, it is expected that it remains almost constant 
since its infall into the host halo. The projected two- 
point correlation fu nction of galaxies is reproduced under 
these assumptions (jConrov. Wechsler fc Kravtsovl 120061 : 
iMoster et al.l l2010h . It should also be said that while 
the (local) SHMR for central galaxies has been deter - 
mined (e.g.. iMandelbaum et alj|2006t iMore et al.ll201ll ), 
the stellar-subhalo mass relation for satellites/subhalos, 
SSMR, has not been yet discussed in detail in the liter- 
ature. 

In view of the above mentioned, some important ques- 
tions arise. Why does using m§££ instead of lead 
to the correct clustering of galaxies? Does the to^-M* 
relation reproduce the observed satellite GSMF, the con- 
ditional stellar mass function, and spatial clustering of 
galaxies at the same time? Even more fundamentally, if 
is not assumed that the SSMR is identical to the central 
SHMR, then, what follows for the SSMR, either using 
or m|£b? Does it deviate from the central SHMR? 

In this paper we extend the common AMT to con- 
strain both the central SHMR and the SSMR separately, 
as well as the average (total) SHMR. By construction, 
this formalism also allows to predict the mean satellite 
conditional mass function (CSMF), i.e., the probability 
that satellites of a given stellar mass reside in distinct 
host halos of a given mass. We will (i) test whether the 
SSMR and the central SHMR follow the same shape; (ii) 
discuss the consequences of defining the subhalo mass at 
accretion time vs. at observed (present) time; and (iii) 
check the self-consistency of our predicted present-day 
central SHMR and SSMR by comparison with the ob- 
served satellite CSMF and the spatial clustering of galax- 
ies. 

This paper is laid out as follows. In Section 2 we 
present the AMT, focusing on the details of our ex- 
tended abundance matching. In Section 3 we present the 
predicted stellar-halo mass relations (§§3.1) and satel- 
lite CSMFs (§§3.2) for cases when the SSMR is assumed 
equal to the central SHMR, and when both mass re- 
lations are independently constrained. In §§3.3, a Halo 
Occupation Distribution (HOD) model is used to explore 
whether the predicted central SHMR's and SSMR's are 
consistent with the observed spatial clustering of galax- 
ies. Section 4 is devoted to our conclusions and a discus- 
sion of the results and their implications. 

All our calculations are based on a flat ACDM cosmol- 
ogy with fl A = 0.73, h = 0.7, and a$ = 0.84. 

2. THE ABUNDANCE MATCHING TECHNIQUE 



1 In order to make the distinction explicit, we shall use upper- 
case letters for the central galaxy and the distinct halo masses and 
lower-case letters for the satellite galaxy and subhalo masses. 



In this section we describe the technique of matching 
abundances between central galaxies and halos and satel- 
lite galaxies and subhalos, separetely, which we present 
here as an extension to the standard AMT. 

2.1. Modeling the central & satellite GSMFs 

To model the central GSMF, let P con (M*|M h ) denotes 
the probability distribution function that a distinct halo 
of mass Mh hosts a central galaxy of stellar mass M*. 
Then the number density of central galaxies with stellar 
masses between M* and M* + dM* is given by 

/>oo 

</> ccn (M*)dM* = dM* / P ccn (M t \M h )(j> h (M h )dM h . 



For the population of satellite galaxies in individual 
subhalos, let P S at("i* |w su b) be the probability distribu- 
tion function that a subhalo m su J3 hosts a satellite galaxy 
of stellar mass m*. Thus the average satellite CSMF (the 
number of satellite galaxies of stellar mass between to* 
and to* + dm,* that reside in distinct host ha los of mass 
M h , e.g. JYang. Mo. fc van den Bosch|[2009bl) is 

$ s (m*\M h )dm* = dm* / P aa t{m*\m snh )$ anh (m snh \M h )dm suh , 



( 2 ) 

where ^sub^subl-Mh) is the subhalo conditional mass 
function (subhCMF, i.e., the number of subhalos of mass 
between m su b a nd TO su b + dm su b residing in host halos 
of mass Mh; e.g.. lBovlan-Kolchin et aLll2010D . A natural 
link between the satellite GSMF, </> sa t, and the distinct 
halo mass function (HMF, 0h) arises once the satellite 
CSMF is given: 

/>oo 

4> sa t{m*)dm* = dm* / ^ s (m*\M h )^ h (M h )dMh. (3) 
Jo 

Inserting equation ^ into equation ^ and rearranging 
terms, the satellite GSMF can be rewritten in terms of 
P sat (m >t |TO su b): 

/>oo 

<p aa ,t(m*)dm* = dm* / P sa ,t(m*\m auh )(f) suh (m suh )dm SVLh , 
Jo 

(4) 

where the subhalo mass function (subHMF) is given by 



4>sub(m sub )dm suh = dm suh / $ sub (TO sub |A/h)0h(Mh)dM h . 

Jo 

(5) 

Equations ((T|) and (U|) describe the abundance matching 
in its differential form for the central-halo and satellite- 
subhalo populations, respectively. The distribution prob- 
ability P cen (M*|Mh) is defined by the mean M*(Mh) re- 
lation and a scatter around it of o~ c , while the distribution 
probability P S a,t{m*\m su b) , assumed to be independent of 
host halo mass, is defined by the mean TO*(m su b) rela- 
tion and a scatter around it of cr s . Observe that once 
Psat(w*|TO su b) is given, the satellite CSMF is a predic- 
tion according to equation @ . 

Here, P CGn (M*|Mh) and Paa.t{m*\ m sub) are modeled 
as lognormal distributions with a width (scatter around 

2 Whenever we use m su ^ we refer to subhalo mass generically. 
In practice, that can either be the mass at accretion time, m-i ?^, 
or at observation (present-day) time, m° i. 



The satellite-subhalo mass relation 
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the stellar mass) assumed to be constant and the same 
for both centrals and satellites, <t c — a s = 0.173 dex. 
Such a value was inferred for central gal axies from the 
analysis of general large group catalogs (1YMB09I ) and 
it is supported b y recent studies o n the kinematics of 
satellite galaxies (jMore etalll20Tl . Re garding the in- 
trinsic scatter of the satellite-subhalo relation, it has 
not been discussed in detail in the literature. While 
the exploration of this scatter is beyond the scope of 
the present paper, our conclusions will not depend criti- 
cally on the assumed value for it or even if it is allowed 
to depend on host halo mass. We will further discuss 
this question in §4.2. Both TO,»(TO su b) and M*(Mh) are 
parametrize d by the same modified two-power-la w form 
proposed in iBehroozi. Conrov fc Wechslerl (|2010T) . This 
five-parameters function is quite general and, in the con- 
text of the AMT, has been shown to reproduce the main 
features of a Schechter-like GSMF. 

2.2. The relation to standard abundance matching 

In the standard AMT the cumulative halo+subhalo 
mass function and the total observed cumulative GSMF 
are matched to determine the mass relation between 
halos and galaxies, which is assumed to be mono- 
tonic. In this context, no intrinsic scatter in the stel- 
lar mass at a given halo is assumed. In our ap- 
proach, where the galaxy and halo populations are sepa- 
rated into centrals/satellites and distinct halo/subhalos, 
the latter entails that the probability distribution func- 
tions of centrals and satellites take the particular forms: 

-Pen 

(M\M h ) = 5(M - M*(M n )) and P sa t(-MKub) = 
5(M — m.*(m su b)), where M*(Mh) and TO*(m su b) are the 
mean central-halo and satellite-subhalo mass relations, 
and 5 is for the (5-Dirac function. The above "no scat- 
ter" probability distribution function for centrals applied 
in Eq. ([1]) would lead us to conclude that the cumulative 
central GSMF, n ccn (> M*), should match the cumula- 
tive distinct halo mass function, rih(> M n (M*)). The 
same reasoning applies for satellites/subhalos. There- 
fore, we arrive to the standard AMT formulation: 



n g (> M*) = n h (> M h ) + n sub (> M h ), 



(G) 



where n g (> M*) = n con (> M*) + n. sa t(> M*) is the total 
GSMF. 

Since the abundance matching can be applied to cen- 
trals/halos and satellites/subhalos separately, let analyze 
now only the latter. Under the assumption that the 
m*(m su b) relation is independent of the host halo mass, 
it is clear that using either the abundance matching of all 
satellites and all subhalos, n sat (> to*) = rt S ub(> TO su b), 
or the matching of their corresponding mean occupa- 
tional numbers, one may find exactly the same to* (TO su b) 
relation. In this sense, we state that matching abun- 
dances is equivalent to matching occupational numbers: 

(N s (> m*|M h )) = (iV BU b(> TO su b | M\i ) ) (7) 
n sat (> to*) = n sub (> m sub ). 

In the case that the probability distribution function 
-Psut^^lTOsub) includes scatter around the mean SSMR, 
as we consider here, the above conclusion remains the 
same whilst P sa t(TO*|TOsub) is assumed to be indepen- 
dent on halo mass. In general, the inclusion of constant 



scatter in the galaxy-halo mass relations is not a con- 
ceptual problem for the AMT, but it slightly modifies 
the shape of the mass relations at the high mass end 
(see IBehroozi. Conrov fe Wechslerl 120101 ) . Finally, note 
that if P sa t(ro* |m su b) depends on Mh, then (j) sa ± may not 
be directly related to </> S ub (see Eq. 4) and using either 
the matching of satellites and subhalo abundances or the 
matching of their corresponding occupational numbers 
would not lead to find exactly the same TO*(m sub ). 

2.3. Inputs for matching abundances 

The inputs required for the procedure described above 
are the subhCMF, the distinct HMF, and the observed 
satellite and central GSMFs. 

For the subhCMF , we use the results obtained in 
iBovlan-Kolchin et al.l (|201dlB~K10D based on the analysis 
of the Millennium-II Simulation. This is a re-simulation 
with the same resolution of a smaller volume of the Mil- 
lennium Simulation. It consists of 2160 3 particles, each 
of mass to = 6.885 x 10 6 /i _1 M Q in a periodic cube of 
length L = 100/i~" 1 Mpc. Observe that this mass parti- 
cle resolution is around four orders of magnitude below 
the subhalo masses required (~ 1 10 /t~ 1 M ( T ) ) to match 
the lower stellar mass limit in the IYMB09I GSMF. The 
fitting formula for the cumulative subhCMF reported in 
IBK10I at the [10 12 , 1O 12 - 5 ]/i _1 M mass interval is: 



{N suh (> m 6 



ib|M h )) = Mo ( ii 



exp 



V Mcut / 



m suh /M h . 



where /i = 
(MOj Ml>Mcut, a, b, 
while for TO sub 

15 (logM h -i2.25) ;0-01) o.096, 

to IBKIOI the shape of the 
the same for other halo masses 
tion, /io, systematically increases 
by 15% per dex in Mh. Such a 



For 



m-sub : 



m 



(8) 

acc 
sub ' 



(1,0.038,0.225,-0.935,0.75), 

— m obs — 

— "'sub) — 

-0.935,1.29). According 



rn 



obs 



sub 



subhCMF remains 
but its normaliza- 
with Mh, roughly 
behavior has been 
reported in an an alysis of the Millennium simulations by 
iGao et all (|2011l for closely related re sults, with small 
differences in the amplitude, see also IGao et all l200l 



van den Bosch. Tormen fc Giocolil 12005 : IZentner et all 
n2009l : IGiocoli etaLl I2010h 



120051: lAngulo et al 

introduce the dependence /io 
where /io = 1 at log Mi, = 
of the m™l subhCMF, the 
nearly independent on Mh, i.e., fin = i i 
IGiocoli. Tormen fc van den Bos ch 2008). 



we 

_ 115 (logJl/ h -12.25) 

12.25. In the case 
normaliz ation factor is 
At n = 1 (|BK10l : see also 
The subhCMF 

is given by $ sub = dN sub /dm s __ 

In order to construct the subhCMF. IBKIOI traced 
each subhalo back in time so that they were able to find 
the point at which its bound mass reached a maximum, 
i.e., the time the halo became a subhalo. The latter 
guarantees that we are working with the surviving pop- 
ulation of accreted halos and no further assumptions on 
the merging process are necessary. 

The difference between the Millennium-II simulation 
cosmology and ours leads to differences in the result- 
ing abundances of subhalos of roughly a few percen t 
in the amplitude of the subhalo mass function (jBKlOh . 
This is also supported by previous works that explored 
the impact of changing cosmological parameters on the 
subhalo occupational statistics (e.g.. IZentner fc Bullock! 
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2003). Additionally, to be consistent with the same cos- 
mology for which the subhalo subhCMF was inferred, we 
repeated all the analysis to be showed below but using 
the WMAP1 cosmology. We find that all our results are 
practically the same. 

F or the distinct HMF, we will use the formula given 
by iSheth fc Tormenl (|1999fl . This formula provides 
a reasonable fit to the the virial mast0 function at 
meas ured in large cosmological N-body sim- 
iKlypin. Truiillo-Gomez fc PrimacH 120111: 



z 

ulations 



iCuesta et al.ll200S 

For our purpose, the decomposition of the GSMF and 
the CSMFs into centrals and satellites gala xies is nec- 
essar y. Using a large general group catalo g (|Yang et all 
[20071) based on the data from the SDSS, YMBOl con- 
structed and studied the decomposition of the GSMF 
and the CSMFs into centrals and satellites galaxies. In 
that paper, a central galaxy was defined as the most 
massive galaxy in a group and the remaining galaxies as 
satellites. For the mass completeness limit in the GSMF, 
they adopted the value as fu nction of redshift proposed in 
Ivan den Bosch etlrtl ([2008). They have also taken into 
account incompleteness in the group mass by considering 
an empirical halo- mass com pleteness limit (for details we 
refer the reader to IYMB09I ). 

Where necessary, halo masses are converted to match 
our definitio n of virial mass and stellar masses are con- 
verted to the lChabrieil ([2003) IMF. In particular, [YMBOl 
defined halo masses a s being 180 times th e background 
density, according to iGiocoli et al.l (|2010L see their ap- 
pendix B) these halos are ~ 11% larger than our defini- 
tion of virial mass. 

2.4. Procedure and uncertainties 

We constrain the parameters of the functions proposed 
to describe the central SHMR and SSMR by means of 
Eqs. (p} and ([3]), and by using the Powell's directions 
set method in m ulti-dimensions for the minimization 
(|Press et a 1.1 1 1992 i. Note that in our analysis the re- 
ported statistical errors in the GSMFs, as well as the 
intrinsic scatter in the mass relations are taken into ac- 
count. However, we will not analyze rigorously here the 
effects of uncertainties on the mass relations as well as 
their e rrors. Instead, we remit the reader to previous 
works dMoster et~al . 2010; B ehroozi. Conrov fc Wechslerl 



[20101 iRodriguez Puebla et al.ll2011l) 



iBehroozi. Conrov fc WechslerT i 20101) studied in detail 
the uncertainties and effects on the average SHMR due 
to different sources of error like those in the observed 
GSMFs, including stellar mass estimates; in the halo 
mass functions; in the uncertainty of the cosmological pa- 
rameters; and in linking galaxies to halos, including the 
intrinsic scatter in this connection. These authors have 
found that the largest uncertainty by far in the SHMR 
is due to systematic shifts in the stellar estimates. The 
second important source of uncertainty is due to the in- 
trinsic scatter, that we take into account in our analysis. 
Other statistical and sample variance errors have negli- 
gible effects, at least for local galaxies. According to the 

3 The mass enclosed within the radius at which, according to the 
spherical collapse model, the overdensity of a sphere is A v ; r times 
larger than the matter critical density of the used cosmological 
model; for the cosmology assumed here, A v i r (z = 0)= 97. 



„ -2 - 



i 

x 
<v 
-a 



o 



so 
o 



-4 - 



-5 - 



"•■•in,.. 




central SHMR: 

SSMR=central SHMR: 

Acc. time v 

Obs. time \ 

SSMR: 

Yang et al. 09: 

Central galaxies B 

Satellite galaxies * 
_i i I i i i i I i i i i_ 



9 10 11 

log M. [M ] 



12 



Fig. 1. — Satellite GSMFs calculated under the assumption that 
Psa.t(jn*\m Euh ) = P CC n(M*\M h ) and for the cases mj°£ (long- 
dashed line) and (short-dashed line) were used for the subhalo 
mass defi nition. Filled circles and squares with error bars show the 
YMB09 central and satellite GSMFs, respectively. The sol id line is 
for the case when P sat (m* |m sub ) was determined using the YMB09 
satellite GSMF as a constraint, i.e., is the best model fit to this 
function. 



IBehroozi. Conrov fc Wechslerl ([2010) study, the statisti- 
cal and systematical uncertainties account for la errors 
in the SHMR of approximately 0.25 dex at all masses, 
which is almost totally due to the uncertainty in stellar 
mass estimates. We have explored here also the effects of 
the subhalo CMF uncertainty on the SSMR . By using the 
25% p er dex in Mh variation reported by IGiocoli et all 
(2010) (instead of 15%), we find that the SSMR shifts in 
m* by only « 0.04 dex. 

3. RESULTS 
3.1. The satellite GSMF and the SSMR 

By means of the procedure described in Section 2, 
we calculate first the satellite GSMF (Fig. [1]) when the 
SSMR and the central SHMR are assumed to be the 
same, i.e., m*(m su b) = Af*(Mh). This is equivalent to as- 
sume that P S at(Ti* |m su b) = P C en(-^* l-Mh) if the intrinsic 
scatter of both relations is the same. We obtain the cen- 
tral SHMR by matching abundances of YMB09 central 
galaxies to distinct halos. This relation and the subhalo 
mass function obtained from the theoretical subhCMF 
(eq. [8]), are used to infer the satellite GSMF (eq. EJ). The 
satellite GSMF is presented for the two cases of subhalo 
mass definition: S at,acc when m^ b is used (long-dashed 
line), and (/> sa t,obs when is used (short-dashed line). 
The observational results of IYMB0 9I are plotted as well. 

Under the assumption that r7i*(m su b) = M*(.Mh), the 
predicted number density of satellites at masses below 
the knee is underestimated on the average by a factor of 
~ 2 when using m^ b , and ~ 5 when u sing m°^. Note 
that the former is closer to the IYMB09I data. The reason 
is simply because the normalization of the fn?££ subHMF 
is higher and closer to the distinct HMF than the nor- 
malization of the subHMF. Therefore, satellites of 
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Fig. 2. — Left panel: Inferred mass relations for satellite galaxies when the subhalo mass is defined as "2°^ ( Dme solid line) and as mJ°E 
(red long-dashed-line), and for central galaxies (black dot- long-dashed line). The density- weighted average relation when using subhalo 
mass w°^S is plotted with a blue short-dashed line, while when using subhalo mass mf* c E is plotted with a gray short-dotted line. For 
comparison, the (average) mass relation obtained in Bchroozi, Conroy & Wcchsler (2010) is also plotted (filled circles). Right panel: A 
comparison of the mass relations of satellite and central galaxies with direct observational inferences (the same line code of left panel for 
models is used; dashed lines indicate extrapolations to lo wer masses). Filled circles with error bars correspond to the mass relation of 
central galaxies from the analysis of staked weak-lensing in Mandclbaum ct al. (2006). Dashed area indicates the 68% of confidence in the 
mass relation of central galaxies us ing the kinematics o f sat ellites lIMore et al.M 2011l) . The inferred total mass at the tidal radii for the 
brightest dwarf galaxies obtained in Strigari ct al. (2007) and lStrigari, Frenk fc White! 1 12010) ar e plotted with fi lled squares and triangles, 
respectively. Filled pentagon shows the mass at the tidal radius for the Large-Magallanic Cloud (Weinberg 1999). The gray dotted-dashed 
curve is the SSMR assuming a faint-end slope in the satellite GSMF of a = —1.6. 



equal to* are expected to have a higher number density 
when using the accreted-time (to^) subHMF compared 
to using the observed-time (mj; present-day) subHMF. 

However, neither mj^£ nor to°^ are able to repro- 
duce the observed satellite GSMF, and the discrepancy 
is due to the basic assumption of a common stellar mass 
(sub)halo mass relation for centrals and satellites. In the 
case that m su b = this is equivalent to assume that 

the SSMR is independent of redshift. But in fact this 
cannot be the case since the satellite mass to* hardly 
will remain the same since it was accreted to the present 
epoch. On the other hand, when using m su b = m°ub> 
that the SSMR is equal to the central SHMR implies that 
both have evolved, on average, identically. This cannot 
be the case because it is evident that the population of 
subhalos evolved differently to distinct h alos, mainly by 
losing mass due to tidal striping (e.g.. iKraytsov et al.l 
[200llvan~den Bosch. Tormen fc Giocolill2005t) . 

The next step in our analysis is to allow the SSMR and 
central SHMR to be different, i.e., TO*(TO su b) ^ M*(Mh). 
In this case, -P S a. t(m*|TO s .i b) is determined by means of 
Eq. ((5) using the lYMB09T satellite GSMF as a constraint. 
For the subhCMF, we again use both definitions of sub- 
halo mass, TOg^ and to°^. For illustrative purpose, we 
present the resulting satellite GSMF for the case when 
m °ub was used (solid line in Fig [TJ an almost identical 
GSMF is obtained when to^ was used). 

As shown in Fig. [21 the SSMRs obtained by using 
(long-dashed line) and to°^ (solid line) are quite differ- 
ent. The central SHMR (dot-dashed line) is the same for 
both cases. The error bar in the left panel shows a lcr 



uncertainty of 0.25 dex in the normalization of the mass 
relations. This is roughly the uncer tainty estimated by 
iBehroozi. Conrov k, Wechslerl (|2010f ) taking into account 
all the systematical and statistical sources of errors (see 
§§2.4). " 

When using the accretion-time subhalo mass, to^,, we 
note that the resulting SSMR at log(m t /M0) < 11 sys- 
tematically lies above the central SHMR, with differences 
in the stellar-mass axis (halo-mass axis) of ~ 0.5 dex (0.2 
dex) at the smallest masses. For log(M*/M©) > 11 this 
trend is inverted, but the differences between central and 
satellites are very small. However, the relation obtained 
this way should be taken with caution. By construction, 
each TOg^ is itself a cumulative distribution of all the 
objects accreted in a time interval Az. Therefore such 
a SSMR entails that all accreted objects of mass mjj££ 
would evolve, on average, to host the same to* despite 
having been accreted at different times. We discuss this 
in §§4.1. 

When using the observation-time (present-day) sub- 
halo mass, to°^, the SSMR (solid blue line) and cen- 
tral SHMR are very different, though they show the 
same trend as when using For example, on av- 

erage, a satellite with log(A/*/M Q ) = 10 resides in a 
subhalo a factor of ~ 4 less massive than the halo of a 
central galaxy with the same stellar mass. Notice that 
m °ub( m *) < m sub( m *) ail d that the difference increases 
the lower the mass is. This is consistent with the picture 
that most massive subhalos, on average, fell into larger 
halos just very recently and they have not had time to 
lose significant amounts of mass due to tidal striping, in 
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contrast to the lowest mass subhalos. 

This also suggests that the SSMR for both definitions 
of subhalo mass should tend to the central SHMR at 
the high-mass end, but this is not the case as seen in 
Fig. [5] where small differences remain. The possi- 
ble reasons are that, firstly, the intrinsic scatter around 
the stellar-(sub)halo mass relations is actually lower for 
the former than for the latter (here we assumed it to 
be the same fo r satellit es and centrals, see §4.2). Sec- 
ondly, that the IYMB09I satellite GSMF may undercsti- 
mat e the true satellite m ass function at large masses (see 
also ISkibba et alJ l20llh . 

Fiber collisions co uld introd uce some systematic error 
that may affect the YMB09| group catalog . To stu dy 
the impact of this possible systematic error, YMBjH di- 
vided their group catalog into two sample: one that uses 
galaxies with known redshifts, and another that includes 
galaxies that lack redshifts due to fiber collisions. When 
compared the corresponding satellite CSMFs from both 
samples (see their Fig. 6), they found that the sample for 
which the correction for fiber collisions has been taking 
into account, has a higher amplitude of the CSMFs than 
when this correction has not been applied, specially in 
low mass halos. However, the difference is very marginal 
and well within the erro r bars. We conclude that fiber 
collisions in the IYMB09I group catalog are not a serious 
source of systematics able to affect our conclusions. Re- 
garding completeness a nd contamination of their group 
catalog (for details see lYang et all 120071 ) , 80% have a 
completeness greater than 0.6, while 85% have a con- 
tamination lower than 0.5. In terms of purity, their halo- 
based group finder is consistent with the ideal situation. 

Finally, we note that the mass relation usually ob- 
tained by matching abundances between the total GSMF 
and the halo plus subhalo mass function, in the light of 
the decomposition into centrals and satellites, could be 
interpreted as a density-weighted average SHMR: 

(M,(M)U = ^^m,(M) + 4^M,(M), (9) 



4>bm(M) 



4>dm(M)' 



where </)um{M) = sub (M) + h (M), m*(Af) is the mean 
SSMR and M*(M) is the mean central SHMR. This re- 
lation is plotted in Fig. [2] with short-dashed-dot and 
short-dashed lines when using m^b an d TO °ub' respec- 
tively. Since most galaxies in the YMB09 catalog are 
centrals, the central SHMR is very close to the density- 
wei ghted average SHMR. For comparis on, we plotted 
the iBehroozi. Conrov & Wechslerl (|2010l ) average mass 
relation (filled circles), which is in excellent agreement 
with our density- weighted average SHMR when using the 
accreted-time subhalo mass, m^. 

Observe that differences between the satellite and the 
average (total) mass relations are small when is 
used, while differences become dramatic when is 
used. The above explains why under the assumption 
that m*(m^ b ) = M*(M h ) = (M*(M)) , the resulting 
satellite GSMF are closer to observations. On the other 
hand, since the m?£Z subHMF has a higher normaliza- 
tion than the subHMF, the above shows that when 
assuming m.fmj) = M*(M h ) = (M*(M))$, we should 
expect that the resulting satellite GSMF is significantly 
below the observed satellite GSMF. 



3.1.1. Comparison with other observational inferences 

In the right panel of Fig. [2 we plot some observa- 
tional inferences of halo and subhalo masses as a func- 
tion of stellar mass. The inferred (Mh)(M*) of central 
galaxies from staked weak -lensing studies using the SDSS 
(|Mandelbaum et"aL 2006) are s hown as filled circles with 
error bars. iMa ndelbaurn et al.l (|2006D reported the data 
actually for blue and red galaxies separately. We es- 
timated the average mass relation for central galaxies 
as: (M h )(M*) = / 6 (M*)<M h ) b (M*) + / r (M*)<M h ) r (M*), 
where fb{M*) and / r (M*) are the blue and red galaxy 
fractions in the sample, and (Mh)& and {Mh) r are the 
corresponding blue and red mass relations. The inferred 
(log(Af*))(Mh) fo r central galaxies from staked kinemat- 
ics of satellites ([More etaJ] l20H are plotted as the 
dashed area indicating the 68% of confidence. 

Our inferred central SHMR (dotted-dashed curve) is 
consistent with the weak-lensing inferences at all masses, 
and with the satellite kinematics inferences at masses 
M* ^ 10 11 M Q ; for smaller masses, our halo masses 
are a factor up to ~ 2 smaller than the satellite kine- 
matics inferences. In fact, it was already noted that 
using the kinematics of satellite galaxies yields halo 
masses around low mass galaxies that are systemati- 
cally larger than most other methods, specially for red 
central galaxies dMore et al.l [20TTt ISkibba et all 120111: 
iRodriguez-Puebla et al.ll2011l ). 

Regarding satellites, unfortunately, there are not di- 
rect inferences of their subhalo masses. Some model- 
dependent estimates based on dynamical observations of 
Milky- Way (MW) satellites were presented in the litera- 
ture. For example, using the line-of-sight velocity disper- 
sion s measured for the br ightes t spheroidal dwarf galax- 
ies. iStrieari et alJ (|2007t ) and iStrigari. Frenk fc Whitel 
(2010) determined their masses within their tidal radii. 
These dynamical masses, plotted in Fig. [5] (filled squares 
and triangles, respectively), are expected to be of the or- 
der of m°„b- We also plot an estimate of the mass at the 
tidal r adius for the La rge-Magallanic Cloud (filled pen- 
tagon, |\^inber3[T999|). The SSMRs constrained here do 
not extend to the small masses of MW satellites but we 
plot their extrapolations to these masses (dashed curves). 
The gray dotted-dashed curve will be discussed in §§4.3 

3.2. The satellite CSMF 

From the approach described in Section 2, another sta- 
tistical quantity that deserves to be subject of study is 
the satellite CSMF (Eq. O . We calculate the mean halo- 
density-weighted CSMF at the [M hl ,M h2 ] bin as: 



•C; 2 MTO*|M h )0 h (M h )dM h 
fi* 4>h(M h )dM h 



(10) 



This qu antitiy has been inferred from observations by 
IYMB09I . again using their SDSS galaxy catalog (filled 
circles with error bars in Fig. [3]). 

First, we consider again the case assuming m*(m su b) = 
M*(Mh). When m sub is defined at the observatio n time, 
the resulting CSMFs are lower than the |YMB09| CSMFs 
by a factor of ~ 5 in the power-law regime (rou ghly the 
same factor by which <^ sa t,obs is lower than the [YMB09 
observed satellite GSMF). Similarly, when m su b is de- 
fined at the accretion time, the predicted CSMFs in 
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the power-law regime are below the |YMB09| CSMFs by 
nearly the same factor, ~ 2, that </> sa t,acc lies below the 
satellite GSMF. The normalization of the CSMF in- 
creases faster with Mh when is used instead of 
m^b- This is because we allow the subhCMF 
normalization to vary with host halo mass, while the 
m sub subhCMF normalization is independent of host halo 
mass. 

The black continuous (mi) and blue long-dashed 
(m^) lines in Fig. [3] (almost indistinguishable one from 
other) are the predictions when P sa .t(jn* \m su b) has been 
constrained by means of the observed satellite GSMF. 
The agreement of the predict ed satelli te CSMF's, for 
both and w ith the IYMB09I CSMF's is now 

remarkable at all halo mass bins for low/intermedium 
stellar masses. Although, as above, the normalization 
of the CSMF's increases faster when m su b = rn°^ than 
when m su b = rn^, the differences between both cases 
at any mass are less than 0.05 dex, within the error bars 
of the observational data. 

Despite the overall agreement, for halo mass bins lower 
than ~ 10 13 Mq, the number of massive satellite galax- 
ies is overestimated, specially at the lowest Mh bins. A 
possible reason for this is the assumption that the scat- 
ter in P sa t(m* |m su b) is constant while in reality it could 
depend on Mh as well as on m su b. However, the prob- 
ability of finding massive satellite galaxies in halos less 
massive than ~ 10 13 M Q is low and they do not con- 
tribute significantly to the mean total density of satellite 
galaxies. Therefore, this assumption does not change our 
conclusions, see also §4. 

Our analysis shows that assuming P sa t(w* m su b) = 
P cen (M Ht |Mh) the resulting satellite CSMFs are not con- 
sistent with observations. Instead, when P sa t(w*|m su b) 
is independently constrained using the observed satellite 
GSMF, there is a clear agreement, no matter what defi- 
nition of m su b was employed for the subhCMF. 

3.3. Abundance matching and clustering 

It has been noted in the literature that the average 
(total) SHMR obtained with the standard AMT is con- 
siste nt with the observed spatial cluste r ing of galax- 
ies (Conrov. Wechsler fc Kravtsovi 120061 : iMoster et al.l 
|2010|. We will test now whether this is the case for 
the mass relations of central and satellite galaxies ob- 
tained here with our extended AMT. We will compute 
the galaxy projected correlation function by means of a 
HOD model for each of the mass relation obtained in 
§3.1. 



A HOD model is a statistical tool mainly 
used to describe th e clu st ering of galaxie s (e.g. , 



Berlind & Weinberg! 120021: ICoorav & Shethl 


2002; 


Yang. Mo & van den Boschl 120031: iZehavi et al.l 


2005, 


20111: iLeauthaud et all l2011allbl: lYang et al.l 12011 


, and 



more references therein). In contrast to the AMT, 
which is a quasi empirical tool, a HOD employs mod- 
eling motivated by resul ts of cosmological A-body 
(e.g.. iKravtsov et al.l 120041 ) and hydrodynamical (e.g., 
Zhe ng et al.ll2005D simulations. 

In short, a HOD model describes the probability that a 
halo of mass Mh hosts a number of N galaxies with stel- 
lar masses greater than M* . Once the occupational num- 
bers are defined, the two-point correlation function can 



be computed assuming that the total mean number of 
galaxy pairs is the contribution of all pairs coming from 
galaxies in the same halo (one-halo term) and pairs from 
different halos (two- halo term). For a detailed descrip- 
tion for the HOD model we employ here, see Appendix 
A. 

First, consider the case when P sa t(w* |m su b) = 
P cen (M*|Mh). The short-dashed curves in Fig. [4] show 
the projected correlation functions in five stellar mass 
bins for the case the subhCMF was used. The 

lYang et"aLl (|2011[) galaxy projected correlation functions 
from the DR7 SDSS are plotted as filled circles with er- 
ror bars. The resulting correlation functions are clearly 
below observations, mainly in the one-halo term. This is 
because using underestimates the satellite GSMF 
and CSMF, resulting in a strong deficit of satellite galax- 
ies. Observe that if N s ~ 0, then N ~ N c and therefore, 
bg(M^) ~ (6(Afh)) where ( fr(M h)) is the mean weighted 
halo bias function, see Eq. IA5I 

When using the subhCMF for m^ b instead of 
(dot-short-dashed curves), the agreement with the ob- 
served correlation functions is better, though at scales 
where the one-halo term dominates, the predictions are 
still below observations. This is, again, because the satel- 
lite GSMF and CSMF are underestimated in this case. 
We remark that using the (average or total) SHMR ob- 
tained with the standard AMT in the HOD model by 
matching the total GSMF to the total halo+subhalo 
mass function (in the case of ?n^ b ), leads to excellent 
agreement with the observed correlation functions, a re- 
sult that is well known. However in this case the SSMR 
is not constrained, instead it is implicitly assumed to be 
equal to the central SHMR (for m^). With our ex- 
tended AMT, we can explicitly separate both mass rela- 
tions. When they are assumed to be equal and the cen- 
tral SHMR is constrained with the central GSMF, then 
we obtain the predictions already shown, in particular 
the correlation functions. The fact that the predicted 
correlation functions, when m^™^ is used, are close to 
those predicted in the standard AMT (and to the ob- 
served ones) is because the central and average SHMR 
are indeed close, as we discussed in §§3.1, see Fig. [2j 

Thus, under the assumption that P sa t(TO*|m su b) = 
P C en(M*|Mh), the observed clustering of galaxies is bet- 
ter reproduced when the subhalo mass in abundance 
matching is defined as rri^Z rather than rn°^. Never- 
theless, even in the former case, the agreement with ob- 
servations is only marginal. 

We now turn the analysis to the cases where the SSMR 
is not assumed to be equal to the central SHMR. The 
black solid and blue long-dashed lines in Fig. [4] show the 
predicted correlation functions in the cases where either 
^sub or TTig^b were used. Both cases lead to very similar 
results and agree very well with observations. 

Therefore, the HOD model combined with the central 
and satellite mass relations independently constrained 
with the extended AMT, is able to reproduce the observed 
correlation functions, no matter if or are 

used. This successful prediction is a consequence of the 
good agreement obtained between our predicted satellite 
CSMFs and those inferred from observations (§§3.2 and 
Fig. 3). 

4. SUMMARY AND DISCUSSION 
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Fig. 3. — Density— weighted average satellite CSMF in eight halo mass intervals. Red short-dashed and green short-dashed-dot lines are 
for the cases when the central SHMR and the SSMR were assumed to be equal and and rnSjSi were used, respectively. The black 

solid and blue long-dashed lines are again for and 

^sub ' respectively, but in the case the central and satellite mass functions were 
independently constrained by means of our extended AMT (they overlap most of time). Filled circles with error bars show the CSMFs 
inferred from observations by YMB09. Note that their halo masses were converted to match our virial definition. 



3 I I I I I I I I I I I I I I I I I I I I r^=r 

7T 3 ± 9<log[M./(h-*M )]<9.5 „ 9.5<log[M./(h-*M o )]<10 „ 



10<log[Kl./(h- 2 M o )]<10.5 



-1 1 I I I I I I I I I J I II "1 I I I I I I I I I I I I 1 I 
<—> : 10.5<logtM./(h- ii Mj]<ll : f ll<log[M./(h- 2 M )]J:11.5 






Yang+11 
P =P 

sub.acc cen 

P =P 
sub.oba cen 



- 



3 AMT P., 
AMT P., 



_L 



10 1-10 1 

log r [Mpc/h] log r [Mpc/h] 



Fig. 4. — Projected two-point correlation functions of galaxies in five stellar mass intervals obtained with the HOD model using different 
galaxy-halo mass relations obtained with our AMT. Gray short-dashed and green short-dashed-dot lines are for the cases when the central 
SHMR and the SSMR were assumed to be equal and m*? and th|5E were used, respectively. The black solid and blue long-dashed lines 
are again for and m^ b , respectively, but in the case the central and satellite mass functions were independently constrained (they 

overlap most of time). The observed projected correlation functions reported in Yang et al. (2011) are shown by filled circles with error 
bars. 
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In this paper, we extend the AMT in order to con- 
strain both the central stellar-halo and the satellite- 
subhalo mass relations separately, using as an input (i) 
the distinct halo and subhalo mass functions, and (ii) 
the observed central and satellite GSMFs. Our formal- 
ism, by construction, predicts also the satellite CSMFs 
as a function of host halo mass, and when applied to 
a HOD model, allows to predict the spatial correlation 
functions. We present results for the cases when the 
SSMR is assumed to be equal to the central SHMR, 
-Psat(w*|m su b)=P C en(^*|-^h), and when both mass re- 
lations are constrained independently (i.e., it is not as- 
sumed that P sa ,t(ni*\m su \ J )=Pcen(M*\M\ 1 )) . All our anal- 
ysis is carried out for subhalo masses defined at accretion 
time, TOguij, and at the observed time (present day), m°^. 
The main results and conclusions are as follows: 

• Assuming that the mass relation between satellites 
and subhalos is identical to the mass relation be- 
tween centrals and distinct halos (including their 
intrinsic scatters), P s:i t(rn*\m sub )=P cen (M*\M h ), 
the predicted satellite GSMF, CSMFs and pro- 
jected two-point correlation functions lie below 
those obtained from observations for both defini- 
tions of m su b, though the disagreements are small 
when m S ub=m^jj (Figs. [TJ [3j |4j. We conclude 
that assuming P sat (m*|m su b)=-Pcen(-M*|Afh) leads 
to predictions in disagreement with observations, 
specially when is used. 

• When the SSMR is no longer assumed to be equal 
to the central SHMR and instead is constrained 
by means of the observed satellite GSMF, the pre- 
dicted satellite CSMFs and projected correlation 
functions agree in general with observations, both 
for mi and Tt^ b . However, for halo masses lower 
than ~ 10 13 Mq, the number of very massive (rare) 
satellites is over-predicted. 

• The resulting m anb -m* relations when using either 
TO °ub or TO sub are quite different from each other, 
and in each case are different from the central 
SHMR (Fig. For a given stellar mass, the satel- 
lite subhalo mass is smaller than central halo mass, 
and the mass difference is increasing the lower the 
mass is. These differences are dramatic when 

is used. 

• Our density-weighted average (centrals + satel- 
lites) SHMRs are close to the central SHMR when 
either mj or m^jj, is used (central galaxies dom- 
inate in the YMB09 catalog). Such an average 
SHMR coincides conceptually with the one inferred 
from matching the total (centrals-l- satellites) cu- 
mulative GSMF and the halo + subhalo cumulative 
mass function (standard AMT). 

4.1. On the inference of the SSMR and its implications 
for the average mass relation 

The conclusions listed above can be well understood 
by examining the basic ideas behind the extended AMT, 
as we show in Section 2.2. Essentially, matching abun- 
dances of satellite galaxies to subhalos is equivalent 



to matching their corresponding occupational numbers, 
that is: 

{N s (> m*|M h )) = {N sub (> m sub \M h )). (11) 

The opposite is also true: matching their correspond- 
ing occupational numbers is equivalent to matching their 
abundances. This is an important result because it shows 
that once P S at(w*|m su b) (and P ce n(M*|Mh)) is properly 
constrained, we will obtain the correct conditional mass 
functions and consequently the correct spatial clustering 
for galaxies. 

The above means that there is a unique m sub (m*) re- 
lationship for each definition of m su b, which depends 
solely on (N sub (> m su b|Mh)). Because of this unique- 
ness, it follows that the rn°^(m^) and m^(m^) rela- 
tions should be different, and any incorrect assumption 
on each one of these relations will lead to inconsistencies 
in the conditional mass functions and spatial clustering 
of galaxies, as for example those that we have found 
here when P sa t(w*|m su b) was assumed to be equal to 
P cen (M*|Mh). Under this assumption, when m^g was 
used, the inconsistencies were actually small. This is 
because in this case the "incorrect" assumption for the 
SSMR is actually not too far from the "correct" result ob- 
tained when P sa t(m*|m°^) is independently constrained 
(compare dot-dashed and solid green curves in Fig. [2]), 
contrary to what happens when is used. 

It is important to remark that in the standard AMT, 
only the average SHMR is constrained (using the to- 
tal GSMF), leaving unconstrained the SSMR, something 
that on its own in troduces a large unc ertainty in the av- 
erage SHMR (see iNeistein et a"fll201lD . We have shown 
that such average SHMR is conceptually equal to the 
density-weighted average mass relation obtained here 
from the observationally constrained central SHMR and 
SSMR. Therefore, our resulting average mass relation is 
expected to be less uncertain than previous determina- 
tions. On the other hand, this average mass relation is 
expected to be close to the central SHMR because most 
of the galaxies in the used observational catalog are cen- 
trals. 

We conclude that in order to properly infer the SSMR 
and the central SHMR at the same time, and this way 
reduce the uncertainty in the average SHMR, more obser- 
vational constraints than the total GSMF are necessary. 
The most obvious and direct is the GSMF decomposed 
into central and satellite galaxies, something that was 
provided by YMB09. However, observe that, according 
to eq. (fTTj) . the satellite CSMFs or the clustering of 
galaxies, modulo the observational errors, provide obser- 
vational constraints that lead to similar inferences of the 
SSMR, because of the uniqueness of this relation for a 
given well defined (N su h(> m su b|Mh)) (see above). 

Finally, we note that obtaining the SSMR for the sub- 
halo mass defined at the accretion time introduces un- 
certainties due to our ignorance about evolutionary pro- 
cesses of the stellar mass since accretion . This does 
not happen when the SSMR is obtained for both the 
satellite and subhalo masses defined at the same epoch, 
for instance the present time. When matching abun- 
dances for the TOj^b case, the fact that (1) m^ b is itself 
a cumulative distribution of all objects accreted over a 
period of time, and that (2) m* may have changed be- 
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tween accretion and observation, are not taken into ac- 
count. In other words, it is implicitly assumed that the 
satellite stellar mass stops evolving soon after accretion. 
In reality the situation is actually quite complex in the 
sense that, depending on the accretion time and the or- 
bit of the satellites, the evolution of their stellar masses 
is diverse, with some of them early quenched and others 
actively evolving, perhaps in some cases as the central 
ones of the same mass (see e.g.. lWetzel, Tinker fe Conrovl 
20111 for semi-empirical inferences on such a complex- 
ity of galaxy evolution in groups). This diversity in- 
troduces an intrinsic uncertainty on the results. Such 
an uncertainty might be accounted for the probabil- 
ity distribution functions: P(rn* m». acc , z), which gives 
the probability that a satellite accreted at epoch z 
evolves, on average, to the observed satellite m„, and 
-P(™*,acc|TOsub' z )j which gives the probability that a sub- 
halo m^jj, hosts a galaxy of mass m». acc at the time of 
accreti on. Now, the satellite CSMF (Eq.[p can be writ- 
ten as (jMo. van den Bosch fc Whitdhoiofc 



$ s (TO*|M h ) 



P(m» |m*, acc , 20^( m *,acc|TOsub> Z ) X 

$ (™subl M h, z)dm*, acc dm^ ) dz. 

(12) 



Note that in our analysis in §§3.1, we implicitly assume 
that the stellar mass of satellite galaxies does not change 
once they become satellites, i.e. P(m» |m». acc , z) = 
5(m* — m„. acc , z), and that P(?ri ! , iacc |m^5i l , z) is indepen- 
dent of redshift. Thus, the application of the AMT to 
infer the satellite CSMF and the m*-m su b relation for 
subhalo mass defined at its accretion time formally re- 
quires more observational constraints at higher redshifts. 
This is a problem already faced by previous authors (e.g., 
lYang et alJlMtl . 

The above is not the only way to formally write the 
satellite CSMF; it can be written in a way that instead 
of implying knowledge of the change of to* from accretion 
to observation, implies just knowledge on the change of 
the subhCMF between these two epochs. Let us consider 
the distribution function, P acc (m°^l\m^, z), giving the 
probability that halos accreted at epoch z evolve, on av- 
erage, to the observed (present-day) subhalos m"^, and 
the probability distribution function of these subhalos of 
hosting a galaxy of mass to* , P sa ,t(m* \m°^). In this case, 
the satellite CSMF (Eq. [2]) is written as 

$ s (m t |M h ) = y J J P s ,t(m4m^)P acc (m^\m^ b ,z)x 

$(m- c u l\M h ,z)dm°^dm^ldz, 

(13) 

and therefore the satellite GSMF, Eq. 01 is given by 



su b(r<£, z)dm°^dm™Zdz. 

(14) 



subHMF, we write 



*K5S) = / / Pacc(m^|m-,z)0 sub (m-,z)dm- b ^. 



( 15 ) 

This last equation is the abundance matching of accreted 
subhalos to present-day subhalos. Therefore, 



Since the subHMF would evolve into the 



<^at(m*) = J Pat(m»|m^)0 sub (m° u b §)dm° u t. (16) 

This equation is nothing but abundance matching satel- 
lite galaxies to subhalos at the time they are observed. 
Hence, the reason that the satellite GSMF matches the 
subHMF in a more direct way for subhalo masses de- 
fined at the observation time (eq. [T6j) than at the ac- 
cretion time (eq. Il2j) . is that in the latter case the 
unknown P(m* m*. acc , z) and P(m ttax:c \m^, z) "evolu- 
tionary" functions have to be introduced. However, we 
acknowledge that for the former case, our ignorance on 
the scatter around the SSMR is a also potential source 
of uncertainty. All our calculations are under the as- 
sumption that this scatter is the same as the scatter of 
the central SHMR. In any case, even if these scatters are 
different, note that including scatter affects the stellar-to- 
(sub)halo mass relation only at its high-mass end, where 
on average satellites are expected to be accreted recently, 
hence their SSMR and scatter are yet similar to those of 
centrals/halos. 

4.2. On the intrinsic scatter in the SSMR 

A possible source of systematic errors in our analysis 
is the assumption that the intrinsic scatter around the 
SSMR, a s is constant and equal to the scatter around 
the central SHMR. To probe the impact of this assump- 
tion we repeated all our analysis but this time assuming 
o~ s = 0. When comparing the results using o~ s = to 
those obtained based on a s — 0.173 dex, we find that 
they are consistent with each other, and therefore with 
the satellite CSMFs an d with th e galaxy spatial cluster- 
ing measured from the YMB09 catalog. In more detail, 
we find that the resulting CSMF's reproduce observa- 
tions for o~ s — slightly better than for a s = 0.173 dex, 
especially at the massive end. This is because when the 
intrinsic scatter is not taken into account (a s = 0), the 
sha pe of the SSMRs steepens at the m assive end (see 
also lBehroozi. Conrov fe We chsler 2010). Consequently, 
for a given m*, the subhalo mass is larger, and the abun- 
dances of larger (sub)halos is lower in general than those 
of smaller halos. Therefore, the number density of satel- 
lites at the massive end is lower. However, the projected 
correlation functions remain almost the same because the 
probability of finding a massive satellite galaxy in host 
halos less massive than logMh ~ 13 is very low. They 
do not contribute significantly to the mean total density 
of galaxies. Although better models are needed in order 
to give a realistic form for a s , our main conclusions seem 
to be robust to variations in the adopted value for o~ s . 

4.3. Implications for satellite/ subhalo evolution 

The local SSMR obtained for both definitions of the 
subhalo mass, and m°^, are such that at halos 

masses smaller than 2 — 10 x 10 13 M Q and at a given 
galaxy stellar mass, the corresponding subhalo mass is 
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smaller on average than the halo mass of centrals (Fig. 
[2]). This difference increases the smaller the mass is, and 
much more for the subhalo mass defined at the observed 
time (present-day). In the case of m^^, the differences 
might be because the halo mass at the epoch it became a 
subhalo (accretion time) is smaller than its present-day 
counterpart at a given stellar mass and/or because the 
satellite stellar mass increased faster than the central one 
for a given halo mass. In fact, it is difficult to make any 
inference in this case because the abundance matching 
is between local galaxies and (sub)halos at different past 
epochs. In any case, the fact that the inferred mass rela- 
tions for satellites and centrals when m^j£ is used are not 
too different, suggests that the central galaxy-distinct 
halo mass relation does not change too much with time, 
at least since the epoch at which most of the subhalos 
were accreted. 

When m°„£ is used, both abundances of satellites and 
subhalos are matched at the same epoch, the observation 
(present-day) time. In this case the strong difference 
between the satellite and central mass relations can be 
interpreted mainly as the result of subhalo mass loss due 
to tidal stripping. Besides, the smaller the subhalo, the 
larger is the mass loss on average. Probably, the different 
evolution in stellar mass between central and satellite 
galaxies could also play a role for the differences but not 
as significant a role the one related to halo and subhalo 
mass evolution. 

From Fig. [2] one sees that for a given m* , the m°^- 
to-m^ ratio is 0.35-0.40 for the smallest masses up to 
m* ~ 2 x 10 11 Mq. At larger masses, this ratio rapidly 
tends to 1. Therefore, the subhalos of satellites galaxies 
less massive than m* ~2x 10 11 M have lost, on aver- 
age (for all host halo massed), 65-60% of their masses 
since they were accreted. It should be noted that this is 
a rough approximation and the evolution of the stellar 
mass since the satellite was accreted should be taken into 
account, see §4.1. This above result shows us that the 
galaxy- (sub)halo connection for satellite galaxies is far 
from direct; present-day satellites of masses m* ~ 7x 10 8 
M and larger have formed in subhalos that at the time 
they were accreted onto galaxy sized halos were on aver- 
age a factor 2.5-3 larger than at the present epoch. This 
has severe implications for studies aimed to constrain the 
A-CDM scenario at the level of subhalo/satellite distri- 
butions. 

For example, it has been discussed that seeding the 
subhalos in simulations of MW-likc halos by using an 
extrapolation to low masses of the stellar-halo mass re- 
lation obtained by means of the AMT, predicts a MW 
dwarf spheroidal (dSph) luminosity function in agree- 
ment with the observed one. However, the circular ve- 
locities at the maximum (or the masses at the infall) 
of the subhalos associated to the dSphs result signifi- 
cantly larger than inferences from observed kinematics 
(jBovlan-Kolchin. Bullock k Kaplinghatll2011[ ). 

In the right panel of Fig. [2] we have plotted the 
extrapolation to low masses of our SSMRs, both for 



4 The dependence of the satellite subhalo mass loss on host halo 
mass will be explored elsewhere. 



subhalo masses defined at the present day (red line) 
and at the infall time (blue line). The observational 
points in the panel are for MW satellites, which sub- 
halo masses were estimated at their truncation radii. 
Thus, if we assume that these masses are roughly equal 
to the present-day subhalo masses in the A-CDM sim- 
ulations, then the simulated subhalo masses, are 
up to w 10 — 30 times larger than those associated to 
dSphs. If the comparison is done with the extrapolation 
of the average (or central) SHMR, then the differences 
increases by a factor of ~ 3 more (see also Fig. 7 in 
iBovlan-Kolchin. Bullock k Kaplinghatll2011[ ). 

Our extrapolated results show that the discrepancy in 
subhalo mass between MW bright dSphs and A-CDM 
simulations is smaller than previously reported but is 
still significant. Note that for the e xtrapolation, we have 
used the same slope of the IYMB09I satellite GSMF at the 
low mass end, a = —1.25 (Fig. Q]). If this slope steepens 
for smaller masses, for example to a value of a = —1.6, 
then our extended AMT predict the SSMR plotted 
as the gray dotted-dashed curve in Fig. [2 which is 
already consistent with the dynamical estimates. 

The GSMF at low masses may be significantly in- 
complete because of missing low-surface brightness 
galaxies. By taking into account the bivariate dis- 
tribution of stellar mass y e rsus surface brightness, 
iBaldrv. Gl azebroo k k Driver! ((2008) have found evi- 
dence for an upturn in the faint-end GSMF slope (a w 
— 1.6) for a subsample of field SDSS galaxies. More re- 
cently, using the G AMA survey, a slo pe of a w —1.47 
has been reported (jBaldrv et al.ll201ll ). Steep faint-end 
slopes have been also found at hi gher redshift s. For 
instance, using the COSMOS field. iDrorv et all (|2009f ) 
have measured slopes of a ~ —1.7 at all redshifts z < 1 
There are also pieces of evidence that the faint-end slope 
of the GSMF (or luminosity function) changes with the 
environment: in clusters of galaxies it steepens signifi- 
cantl y (for a discussion see IBaldrv. Glazebrook k Driverl 
120081 and the references therein). The cluster GSMF 
is actually related to the satellite GSMF, through the 
satellite CSMF. 

We conclude that using a correct AMT for connect- 
ing satellite galaxies to their present-day subhalos and 
assuming a steep faint-end slope in the satellite GSMF 
(a 1.6), the predicted SSMR for the A-CDM cos- 
mogony would be consistent with the dynamics of MW 
satellites. 
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APPENDIX 

THE SPATIAL CLUSTERING OF GALAXIES IN THE HOD MODEL 

Here we review the main ideas used to infer the spatial clustering of galaxies based on a HOD model. We assume 
that the most massive galaxy in terms of stellar mass within a halo of mass is its central galaxy. Consequently 
the remaining galaxies are all satellites. We let them follow the mass density profile of the host halo. We denote the 
cumulative number of central and satellite galaxies with stellar masses greater than M* as N c and N s , respectively. 

The two point correlation function is decomposed into two terms, 

l + Wr) = [l+CM] + [l + ^ g h M], (Al) 

with 1 + £gg( r ) describing galaxy pairs within the same halo (the one-halo term), and 1 + C|g( r ) describing the 
correlation between galaxies occupying different halos (the two-halo term). 
We compute the one-halo term as 

1 + &V) = ^-1^ J™ {N{N - 1)} X(r)MM h )dM h , (A2) 

for pairs separated by a distance r ± dr/2. Here (N(N — l))/2 is the total mean number of galaxy pairs within halos 
Mh following a pair distribution X(r)dr, and a mean number density n g (M*). The contribution to the total mean 
number of galaxy pairs from central-satellite pairs and satellite-satellite pairs is 

(N(N - 1)) = {Nc){Ns)XcAr)dr 
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+ ™^»A„ W *. (A3, 

As commonly assume d in HOD models, the number o f central-satellite pairs follow the normalized mass density halo 
profile, taken to be of iNavarro. Frenk fc White] (jl997) shape. The number of satellite-satellite pairs is then related to 
the normalized density profile convolved with itself. 

Halo profiles are defined in terms of the total halo mass and t he concentration parameter. We use the relation 
between concentration parameter cnfw an d halo mass obtained by IMunoz-Cuartas et al.l (|2011D from fits to iV-body 
simulations. 

Based on results of high- resolution TV-body ()Kravtsov et al.ll2004f) and hydrodynamic simulations of galaxy formation 
(|Zheng et al.l 120051 ). we model the second moment of satellite galaxies, {N S (N S — 1)}, as a Poisson distribution with 

mean (N s ) 2 = {N S (N S - 1)}. 
We compute the two-halo term as 

&) = blC 2 (r)U(r), (A4) 

where £ m (r) is t he non-linear matte r correlation function (|Smith et al.ll2003|) . £(r) is the scale dependence of dark 
matter halo bias (jTinker et al.ll2005h ■ and 

b g = — [ b(M h )(N(> M*\M h ))MM h )dM h , (A5) 

n 9 Jo 

is the galaxy bias with b(Mh) the halo bias function (Sh eth. Mo fc Tormenl[200ll ). 

Once we have calculated £ gg (r), we relate it to the projected correlation function, u> p (r p ), integration over the line 
of sight, 

w p (r p ) = 2 r ^rl + x^dx. (A6) 



